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Abstract 

We study the nonequilibrium Josephson current in a long two-dimensional ballistic 
SNS-junction with a normal reservoir coupled to the normal part of the junction. 
The current for a given superconducting phase difference (j) oscillates as a function 
of voltage applied between the normal reservoir and the SNS-junction. The period 
of the oscillations is Trhvp/L, with L the length of the junction, and the amplitude 
of the oscillations decays as for eV S> hvp/L and zero temperature. The 

critical current Ic shows a similar oscillating, decaying behavior as a function of 
voltage, changing sign every oscillation. Normal specular or diffusive scattering at 
the NS-interfaces does not qualitatively change the picture. 



pacs[74.50.+r, 74.20.Fg, 74.80.Fp] 



1 Introduction 



In recent years there has been an increased interest in the nonequilibrium 
Josephson current in mesoscopic multiterminal SNS-junctions. Quasiparticle 
injection from one or several normal reservoirs coupled to the SNS-junction 
leads to a nonequilibrium population of the current carrying Andreev levels, 
and thus to a modification of the Josephson current. As predicted by theory, 
[1-6] experiments show that suppression [7,8], switching [9] and even enhance- 
ment [10] of the Josephson current under nonequilibrium is possible. 

The theory for the nonequilibrium Josephson effect has been developed for 
mainly two types of junctions, quantum ballistic [1,3,4] and diffusive [2,6,5]. 
A growing experimental interest is however shown for multiterminal SNS- 
junctions where the normal part is a ballistic semiconductor 2DEG.[8,11] In 
this paper, a theory for these type of structures is presented, and it is applied 



Preprint submitted to Elsevier Preprint 1 February 2008 



to junctions with and without normal reflection at the NS interfaces. For 
junctions with normal reflection, both specular and diffusive scattering is taken 
into account. 



2 Model and theory 

A model of the junction is presented in Fig. 1. A ballistic two-dimensional 
normal region, width W, is connected to two superconducting electrodes, with 
electrode separation L. The phase difference between the superconductors is 

</). A normal electron reservoir is connected to the normal part of the junction 
via a quantum point contact, with width d <^ W, and a voltage V is applied 
between the normal reservoir and the SNS-junction. 



The resistance of the point contact is assumed to be the dominating resis- 
tance of the junction, and the applied voltage thus drops completely over the 
injection point. Zero magnetic field is assumed. 

A natural framework for studying multiterminal ballsitic two-dimensional SNS- 
junctions is the Landauer-Biittiker scattering approach. The junction is de- 
scribed by the Bogoliubov-de Gennes-equation (BdG), where it is assumed 
that the superconducting gap A is constant in the superconductors and zero 
in the normal metal. Hard wall boundary conditions for the 2DEG are also 
assumed. Solutions to the BdG-equation are matched at the NS-interfaces and 
at the three lead connection. It is assumed that all transverse modes couple 
equally to the modes in the injection point contact. [12] 

Boundary conditions are incoming electron and hole quasiparticles from the 
normal reservoir and incoming electron- and hole-like quasiparticles from the 





Fig. 1. A schematic picture of the junction. 
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superconductors at energies above the gap. Knowing the wave function coef- 
ficients the current density is straightforwardly calculated. 

Since the normal reservoir is weakly coupled to the SNS-junction {d <S W), the 
injected current is small and the current flowing between the superconductors 
is only the nonequilibrium Josephson current. [13] However, the coupling must 
be large enough so that the injected quasiparticles do not scatter inelastically 
before leaving the junction. Under these assumptions, the distribution of the 
injected electrons and holes are governed by the normal reservoir. 

The total nonequilibrium Josephson current is naturally parted into the equi- 
librium current /^^ (at eV = 0) and the current due to nonequilibrium, 
It has been shown [14] that the nonequilibrium current can be split into two 
components, one associated with the nonequilibrium population of the An- 
dreev states and the other with nonequilibrium mesoscopic fluctuations of the 
current. Here, this mesoscopic fluctuation term is neglected, since it is small 
compared to term from the nonequilibrium population of the Andreev levels. 
The total nonequilibrium current then becomes 



/ = j^i _|_ j"-^<i 



J dE i{E) nF+ J dE 



2nF) 



with n^*^^) = np{E =]= eV) being the distribution functions of electrons (holes) 
in the normal reservoir, where np = [l-|-exp(£^/A;T)]~^. Clearly, the properties 
of the current density i{E) directly determines the nonequilibrium Josephson 
current, and the current density will thus be the staring point for the discus- 
sions below. For energies outside the gap, E > A, i{E) is given by [15,16], 



4e d 
h d(t) 



Im (tr In [l - a'^S{E)rAS*{-E)r*^) , 



(2) 



where S{E)[S*{—Ey\ is the electron(hole) scattering matrix for the normal 
region and a — exp[— acosh(ii^/A)] and Ta — diag[exp(i0/2), exp(— i0/2)] de- 
scribe the Andreev reflection at the NS-interfaces. The dimensions of the scat- 
tering matrices S{E) and are 2A'" x 2A'", with = 2W/\f the number of 
transverse transport modes in the normal region between the superconduc- 
tors, with \f the Fermi wavelength. The current density for energies within 
the gap is given by adding a small, positive imaginary part to the energy 
E ^ E + i0.[17,16] The expression (2) then reduces to the well known result 
[18] 

m = 'i^'-^HE-E„). (3) 
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where the index m is labehng the bound states given by the equation det[l — 
a"^ S {E)rAS* E)r\] = 0. This form of the current density is useful when the 
bound state energy as a function of phase difference (f) is explicitly known. The 
current density has the energy symmetry i{^E) = —i[E). 

The junctions studied are in the long limit, L ^ ,^oi Co = hvp/^. The 
nonequilibrium Josephson current injunctions in the opposite, short limit, has 
been studied in Ref.[14]. Only classically wide junctions, with many transport 
modes N ^ 1, are considered below. The opposite quantum limit, = 1, was 
studied in Refs. [1,3,4]. 



3 Perfectly transmitting NS-interfaces 



We first consider the case when there is perfect Andreev reflection at the NS- 
interfaces, i.e no normal reflection. For the low energy levels, £J ^ A, the 
bound state energies are given by [19] 

i^;n=^[(2p+lK±0], (4) 



where the index n denotes the transverse mode and p, ± labels the Andreev 
levels for a given mode. Due to the hard wall conditions, the Fermi velocity for 
each mode n\s vpn — vp^j'^- — {n/Ny. The current density for each Andreev 
level is given by inserting the expression (4) into (3). The total current density 
is then obtained by first summing over the modes n, equivalent to integrating 
over angles, and then summing over p, ±, giving 



where 9 is the Heavyside stepfunction and E^ is given from Eq. (4) with 
n = 0. The current density for phase differences (f) — n/A and 37r/4 is plotted 
in Fig. 2. The current density consists of alternating positive and negative 
peaks at energies E^q. The peaks arise from the square root singularities of 
the current density in Eq. (5) and the amplitude of the peaks is decreasing for 
increasing energy. 

The expression for the equilibrium part, of the Josephson current in Eq. (1) 
is known. [19-21] The nonequilibrium part, is straightforwardly calculated 
from Eqns. (5) and (1). It is clear from the form of the current density in Eq. 
(5) that the total current / will be an oscillating function of voltage, with 
alternating maxima and minima, at voltages eV — E^ for zero temperature 
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Fig. 2. Upper: The current density as a function of energy. Lower: The total cur- 
rent as a function of voltage for zero temperature. The phase difference = 7r/4 
(dashed)and 37r/4 (solid) and L = lO^o 

(see Fig. 2). The period of oscillation is thus nhvp/L. The amplitude of the 
oscillations, Alp — I{eV — E^) — I{eV — E'q), in the limit eV » hvp/L, 
decays with voltage as 



N 



evp 



( \(P\^VF 

L \ LeV 



3/2 



(6) 



For finite temperatures the amplitude of the voltage oscillations decreases, 
and the oscillations are completely washed out ior kT ^ kvp/L. 
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The equilibrium current /^^ is always positive for phase differences < < tt 
and negative for — tt < < 0. In nonequilibrium, the total current for a 
given phase difference I{4>) changes sign as a function of applied voltage, i.e 
the junction becomes a so called 7r-junction [2,9]. This can bee seen from the 
current to phase relationship I {(f)), shown in Fig. 3 for different voltages. It 
can be noted that for certain voltages, the current to phase relationship has 
several local current minima and maxima. 

The critical current is defined as the maximum possible current for phase 
difference — vr < < tt. To study the 7r-junction behavior, the critical current 
multiplied by the sign of the critical phase difference (6,., sgn(0c)/c, is shown in 
Fig. 3 for different temperatures. The critical current multiplied by the critical 
phase difference sgn(0c)-^c oscillates between positive and negative values as 
a function of voltage, with a period ttHvf/L, i.e it shows a typical 7r-junction 
behavior. The amphtude of the oscillations decreases with increasing voltage. 
It can be noted that for zero temperature, there are, for certain voltages, jumps 
between the branches of positive and negative critical current, i.e the critical 
current Ic never becomes zero. This can be understood from the current-phase 
relationship in Fig. 3, where for certain voltages, the critical phase difference 
0c jumps between positive and negative values, changing the sign of sgn(0c)-^c- 



4 Normal reflection at the NS-interfaces 



Normal scattering at the NS-interfaces is taken into account by introducing 
effective interface barriers with the transmission probability F. The low lying 
bound states energies, £^ -C A, are given by[22] 



4(1 - F)cos(2A;^„L) - F2cos0\1 
(2 - F)2 JJ ■ 

The bound state energies oscillate rapidly as a function of length of the junc- 
tion, due to the term cos(2/cf„L) in Eq. (7). The corresponding rapid oscilla- 
tions of the current density are averaged out when summing over the transverse 
modes. [23] The total current density, summed over n and p, ±, is given by 



E, 



± 

p,n 



Fn 



2L 



2pn ± acos 



. ^ e 2r2sin0 sg^[sm{2EL/hvFj\ 

^ ^ y ^^16(1 - F)2 - [(2 - Vycos{2ELinvFr) + V^coscff 

The current density as a function of energy is plotted for different barrier 
transparencies F in Fig. 4. The current density has a similar shape, with 
alternating positive and negative peaks, as the current density for the junction 
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Fig. 3. Upper: The current phase relationship for different voltages eV = (dotted), 
eV = hvp/L (dashed) and eV = 2?ivp/L (solid), temperature T = 0. Lower: The 
critical current multiplied by the sign of the critical phase difference sgn(^c)/c as a 
function of applied voltage for different temperatures, kT = 0, 0.3, O.S^vf/L, with 
decaying amplitude for increasing temperature. 

without barriers at the NS-interfaces, shown in Fig. 2. The effect of the NS- 
barriers, apart from the overall decreased current density amplitude, is that 
each current density peak is shifted towards lower energies, as is seen in Fig. 4. 

The current as a function of voltage, for a fixed phase difference, thus oscillates 
with the same period nhvp/L as in the case without barriers at the NS- 
interfaces, and with a decreasing amplitude of the oscillations for increasing 
voltage (see Fig. 4). 
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Fig. 4. Upper: The current density as a function of energy. Lower: The total cur- 
rent as a function of appled voltage. The barrier barrier transparencies F = 0.1 
(solid), 0.5 (dashed) and 0.9 (dotted). The phase difference <p = 37r/4 and the 
length L = lO^o 



5 Dirty NS-interfaces 



Usually, there is dirt at the NS-interfaces from the junction processing, leading 
to diffusive scattering at the interfaces. To simulate this diffusive scattering, 
we introduce random scattering matrices 5*1 and 5*2 to model the interfaces. 
The matrices are written in the polar decomposition [24] 



(9) 



where the barrier transmittances Tj are taken to be mode independent and 
equal for both barriers. The unitary matrices Uj, Vj are taken to be indepen- 
dent members of the ensemble of unitary, symmetric matrices (COE)[25]. The 
total electron scattering matrix for the normal region is given by 



S ^r + tSo{l-r'So)-V, So 




(10) 



where r = diag(ri,r2) and similarily for r' and t. The matrix 5*0 is the scat- 
tering matrix for the normal region without NS-barriers, with the diagonal 
matrices P with elements P„ = exp(i[A;i;'„ -|- EL/{hvFn)]), kpn beeing the 
Fermi wave vector. This is inserted in Eq. (2) to give the current density. It 
can be pointed out that for only forward mode mixing scattering (F = 1), the 
mode mixing at each interface is completely reversed by the Andreev reflec- 
tion, giving the same result as in the absence of barriers at the NS-interfaces. 

We are interested in the current density averaged over the random matrices 
Uj, Vj, which is calculated numerically by generating a large number of matri- 
ces. The current density is plotted in Fig. 5 for different barrier transparencies 
F. 
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Fig. 5. The current density as a function of energy for different barrier transparen- 
cies r = 0.1 (solid), 0.5(dashed) and 0.9(dotted). The current density has been 
calculated by generating 1000 random matrices of dimension N = 15. The phase 
difference <p = Sir/A and the length L = lO^o 
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The main result is that the peak-like structure of the current density is not 
substantially changed compared the corresponding specular barrier case in Fig. 
4. From this one can draw the conclusion that the effect of diffusive scattering 
at the NS-interfaces does not modify the nonequilibrium Josephson current 
quahtatively. The result that the current for a fixed phase difference oscillates 
periodically with applied voltage is still valid. Whether rough 2DEG-sidcwalls, 
giving rise to diffusive boundary scattering, have a more profound effect on 
the current, remains to be investigated. It is known that in the case with the 
normal region being a chaotic cavity, the peak-like current density structure 
is completely washed out, and a gap opens up in the spectrum. [16] 



6 Conclusions 

In conclusion, we have studied the nonequilibrium Josephson current in long 
two-dimensional ballistic SNS-junctions weakly coupled to a normal metal 

reservoir. The total current is given by a convolution of a single current density 
i{E) with the quasiparticle distribution functions [See Eq. (1)]. Junctions with 
and without specular normal scattering at the NS-interfaces and also junctions 
with diffusive NS-interfaces are studied. It is found that the current density 
in all cases has a peak-like structure, with alternating signs of the peaks. 
The resulting nonequilibrium Josephson current for a given phase difference 
thus oscillates as a function of applied voltage, with a period nhvp/L, and 
an amplitude decreasing with increasing voltage. This behaviour also carries 
over to the critical current, which changes sign as a function a voltage, i.e the 
junctions displays so called 7r-behavior. 
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